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Abstract — Joint network-channel codes (JNCC) can improve 
the performance of communication in wireless networks, by 
combining, at the physical layer, the channel codes and the 

[ network code as an overall error-correcting code. JNCC is 
increasingly proposed as an alternative to a standard layered 
construction, such as the OSI-model. The main performance 
metrics for JNCCs are scalability to larger networks and error 

] rate. The diversity order is one of the most important parameters 
determining the error rate. The literature on JNCC is growing, 
but a rigorous diversity analysis is lacking, mainly because of the 
many degrees of freedom in wireless networks, which makes it 

\ very hard to prove general statements on the diversity order. In 
this paper, we consider a network with slowly varying fading 
point-to-point links, where all sources also act as relay and 
additional non-source relays may be present. We propose a 

] general structure for JNCCs to be applied in such network. 
In the relay phase, each relay transmits a linear transform 
of a set of source codewords. Our main contributions are the 
proposition of an upper and lower bound on the diversity order, 

] a scalable code design and a new lower bound on the word 
error rate to asses the performance of the network code. The 
lower bound on the diversity order is only valid for JNCCs 
where the relays transform only two source codewords. We then 

[ validate this analysis with an example which compares the JNCC 
performance to that of a standard layered construction. Our 
numerical results suggest that as networks grow, it is difficult to 
perform significantly better than a standard layered construction, 

I both on a fundamental level, expressed by the outage probability, 
as on a practical level, expressed by the word error rate. 



I. Introduction 

Point-to-point communication has revealed many of its se- 
crets. Driven by new applications, research in wireless com- 
munication is now focusing more on the optimization of 
communication in wireless networks. For example, the joint 
operation of multiple network layers can be optimized, denoted 
as cross-layer design 1361 . 1371 . thereby leaving the classical 
layered architectures, such as the seven-layer open systems 
interconnect (OSI) model fT, P- 20]. Another example of 
network optimization is cooperative communication, where 
multiple nodes in the network cooperate to improve their 
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error performance. Cooperation may occur in many forms at 
different layers, e.g. cooperative channel coding at the physical 
layer and network coding at the network layer Network coding 
refers to the case where the intermediate nodes in the network 
are allowed to perform encoding operations over multiple 
received streams from different sources. In a standard layered 
construction, the decoding of the network code is performed at 
the network layer, after the point-to-point transmissions have 
been decoded at the physical layer Channel coding refers to 
the case where nodes perform coding over one point-to-point 
wireless link only. Cooperative channel coding is achieved 
by letting one or more relays transmit redundant bits for one 
source at a time. Usually, channel coding and network coding 
are studied separately (e.g. ||9l, ED, 1271 for cooperative 
channel coding and HI, El, EH), ||34|, li2 for network 
coding). 

Joitit network-channel coding (JNCC) received much of 
attention in the last years. This technique combines both de- 
code and forward l25]| (cooperative communication) and cross- 
layer design by using a network code, which is accessible 
at the physical layer The rationale behind joint network- 
channel coding is to improve the performance, by combining, 
at the physical layer, the channel codes and the network 
code as an overall error-correcting code ITSl . Mostly, the 
two most important performance metrics are {R,Pe), where 
R is the spectral efficiency and Pg is the error rate (bit 
error rate or word error rate). Here, we consider a fixed 
spectral efficiency R, so that the aim is to minimize Pe for a 
given channel quality, expressed by 7, the signal-to-noise ratio 
(SNR). Expressing the asymptotic (for large 7) error rate as 
Pe = where c and d are defined as the coding gain and 
the diversity order, respectively, improving the performance 
refers to maximizing first d and then c (because d has the 
larger impact). 

Standard linear network coding consists of taking linear 
combinations of several source packets, well known as the typ- 
ical XOR operations for binary codes. In general, non-binary 
coefficients are used in the linear combinations. However, 
when the network code is used at the physical layer to decode 
the noisy channel output, this simple technique might yield 
poor error performance. Therefore, powerful network codes 
consisting of taking linear transformations of the incoming 
information packets, have been introduced. We denote this 
methodology as generalized linear network coding (GLNC). 
The well known standard linear network codes, taking linear 
combinations, are a special case of GLNC. Combining GLNC 



2 



with channel coding, is denoted as joint network-channel cod- 
ing (JNCC). The JNCC, which is the overall code comprising 
the channel codes and the network code, can for example be 
an LDPC code or a Turbo code. Of course, while JNCC brings 
more degrees of freedom and opens perspectives for a higher 
coding gain c, it must be verified that important metrics, such 
as the diversity order d and the scalability to large networks, 
are not negatively affected. 

Binary JNCC has already been studied in the literature. 
Pioneering papers ifTTl . ifTSl designed Turbo codes and LDPC 
codes, respectively, for the multiple access relay channeQ 
(MARC) and for the two-way relay channel |fT9l . However, 
the code design was not immediately scalable to general 
large networks and did not contain the required structure to 
achieve full diversity. In HOl, HIl, a full-diversity JNCC for 
the MARC was proposed but it was not extended to large 
networks. The work of Hausl et al. lfT7l - lfT9l was followed by 
the interesting work of Bao et al. 121, presenting a JNCC that 
is scalable to large networks. However, this JNCC was not 
structured to achieve full diversity and has weak points from 
a coding point of view fT2l . A deficiency in the literature, 
for general networks with sources and relays, is the lack of 
a detailed diversity analysis in the case that the sources can 
act as a relay (which is for example the model assumed by 
||2l). The effect of the parameters of the JNCC on the diversity 
order is in general not known, because it is very hard to prove 
general statements on the diversity order, in an environment 
with so many degrees of freedom. This paper is a modest 
attempt to contribute to the solution of this problem. Related 
to this, we mention 1291 . 1301 . where the authors designed a 
JNCC for the case where the sources cannot act as a relay, 
but other nodes play the role of relay to communicate to 
one destination. As the source nodes are excluded to act as a 
relay node in this model, the diversity analysis in l29l . ll30l 
is different from ours. 

In this paper, we consider a JNCC where the network code 
forms an integral part of the overall error-correcting code, that 
is used at the destination to decode the information from the 
sources. The body of this paper consists of two main parts. 
In Sec. |IV] we perform a diversity analysis, leading to an 
upper bound on the diversity order of any linear binary JNCC 
following our system model, and to a lower bound on the 
diversity order for a particular subset of linear binary JNCCs. 
The upper and lower bound depend on the parameters of the 
JNCC and can be used to verify whether a particular JNCC 
has the potential to achieve full diversity on a certain network. 
Secondly, in Sec. |V] a specific JNCC of the LDPC-type is 
proposed that achieves full diversity for a well identified set 
of wireless networks. The scalability of this specific JNCC 
to large networks is discussed. The coding gain c is not 
considered in the body of the paper and the parameters of 
our proposed code may be further optimized by applying 
techniques such as in l29l . to maximize c. To assess the 
performance of the proposed JNCC, we determine the outage 
probability, a well known lower bound of the word error rate, 

'in the MARC, two sources are helped by one relay to communicate to the 
destination. 



in Sec. |VI] We also present a tighter word error rate lower 
bound in Sec. IVI-BI that takes into account the particular 
structure of the JNCC. In Sec. IVIII the numerical results 
corroborate the established theory. We also briefly comment 
on the coding gain achieved by the proposed JNCC and 
conclusions are drawn for different classes of large networks. 

This paper extends the work, pubUshed in IT2l . by also con- 
sidering non-perfect source-relay channels, by considerably 
extending the diversity analysis, by providing an achievability 
proof for the diversity order of the proposed JNCC, by clearly 
indicating the set of wireless networks where the proposed 
JNCC is diversity-optimal, by providing a tighter lower bound 
on the word error rate, and by providing more numerical 
results. 



II. Joint network-channel coding 

We first illustrate joint network-channel coding by means of 
a simple example. Consider two sources orthogonally broad- 
casting a vector of symbols, mapped from the binary vectors 
Si and S2, respectively, to a relay and a destination. This 
channel is denoted as a multiple access relay channel (MARC) 
in the literature. Supposing that the relay is able to decode the 
received symbols, the relay computes a binary vector ri, which 
is mapped to symbols and transmitted to the destination. The 
relation between all bits is expressed by the JNCC, whose 
parity-check matrix has the following general form, 

si S2 ri 

-Hp - 

H= \ I ■ (1) 

H<y Hi 

The matrix Hp represents the parity-check matrix for the point- 
to-point channel code. Each of the binary vectors Si, S2 and ri, 
can be separately decoded using this code. The bottom part 
of H represents the GLNC, which we denote as i?GLNC = 
[Hi H2 Hi]. It expresses the relation between ri, Si and S2. 
More specifically, we have 

Hin = Hlsi + Hls2. (2) 

Note that GLNC includes standard network codes used in an 
OSI communication model as a special case. In the latter case, 
the matrices i/j and Hi (considering more than one relay 
in general) are identity matrices or all-zero matrices, so that 
the network code simplifies to the relay packet being a linear 
combination of source packets, also expressed as XORing of 
packets or symbol-wise addition of packets. 

Ideally, the overall matrix H conforms optimized degree 
distributions that specify the LDPC code. When the channels 
between sources and relay are perfect, we can drop the first 3 
sets of rows and only keep the GLNC, represented by iJcLNc; 
in this case the information bits of the code are Si and S2, 
and ri contains the parity bits. This is still a JNCC as the 
redundancy in the network code is used to decode the received 
symbols on the physical layer at the destination. In ifTol . ITTl . 
it is proved that the matrices Hp do no affect the diversity 
order in the case of the MARC. 
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III. System model 

We consider wireless networks with sources directly 
communicating to a common destination (e.g. cellphones com- 
municating to a base station). Two time-orthogonal phases are 
distinguished. In the source phase, the sources orthogonally 
broadcast their respective source packet. In the following 
relay phase, the relays orthogonally broadcast their respective 
packet. All considered sources overhear each other during the 
source phase, and act as relay in the relay phase. Other nodes, 
not acting as a source, might be present in the network (i.e., 
overhearing the sources) and also act as relay. Hence, we 
consider a total of nir relays, where > rag. This general 
network model, which is practically relevant as it fits many 
applications, is adopted in e.g. 121. Take for example any 
large network and consider a volume in space (cfr. picocells 
or femtocells) where all nodes can overhear each other. These 
nodes form sub-networks and can be modelled by our proposed 
model. Note that in the literature, sometimes other models are 
assumed, such as the M — N — 1 model ||29l , ll30l . where M 
sources are helped by N relays (the relays are nodes different 
from the sources) to communicate to one destination. 

All devices have one antenna, are half-duplex and transmit 
orthogonally using BPSK modulation. The K information bits 
of each source are encoded via point-to-point channel codes 
into a systematic codeword, denoted as source codeword, of 
length L, expressed by the column vector for user Us,Us G 
[1, . . . , rus]- The parity-check matrix of dimension (L—K) x L 
of this point-to-point codeword is denoted by Hp, which is 
the same for each user Ug, so that i?pS„^ = for all Ug. In 
the relay phase, each relay Ur, G [1, . . . ,TOr]; transmits a 
point-to-point codeword r^^ of length L to the destination, also 
satisfying HpYu^ = 0. Hence, all slots have equal duration, the 
coding rate of the point-to-point channels is i?c.p = ^, and 
the overall coding rate is Rc = -, — s , = Rr-o — ^r^ — ■ 
We define the fraction of source transmissions in the total 
number of transmissions as the network coding rate i?„ = 
— — , so that Rn = Rn r,Rn- The overall codeword of length 



{nig + mr)L is expressed by the column vector 



(3) 



The destination declares a word error if it can not perfectly 
retrieve all mgK information bits, and the overall word error 
rate is denoted by Pew- 
All relevant channels between different pairs of network 
nodes are assumed independent, memoryless, with real ad- 
ditive white Gaussian noise and multiplicative real fading 
(Rayleigh distributed with expected squared value equal to 
one). The fading coefficient of a wireless link is only known 
at the receiver side of that link. We consider a slow fading 
environment with a finite coherence time that is longer than 
the duration of the source phase and the relay phase, so that 
the fading gain between two network nodes takes the same 
value during both phases. We denote the fading gain from 
node u to the destination as a^, with E[a,^] = 1. All point- 
to-point channels have the same average signal-to-noise ratio 

^Unless mentioned otherwise, we assume that channels are reciprocal, i.e., 
the channel from ui to U2 is the same as the channel from U2 to ui. 



(SNR), denoted by 7. Differences in average SNR between 
the channels would not alter the diversity analysis, on the 
condition that the large SNR behaviour inherent to a diversity 
analysis refers to aljj SNRs being large. Denoting the received 
symbol vector at the destinatior0 in timeslot i as y^, the 
channel equation is 



(4) 



Yu, = aM,s„^ + n„^, Ug^l,...,mg 

where ^ CA/'(0, -/) is the noise vector in timeslot i, = 
2s„^, - 1 and = 2r„^ - 1 (BPSK modulation). 

Hence, at the destination, each of the irig independent fading 
gains between the sources and the destination affects 2L bits 
(L bits in the source phase and L bits in the relay phase) and 
each of — nig fading gains between the non-source relays 
and the destination affects L bits, assuming that all nir relays 
could decode the messages received from the sources. Hence, 
from the point of view of the destination, the overall codeword 
is transmitted on a block fading (BF) channel with blocks, 
each affected by its own fading gain, where irig blocks have 
length 2L and m^ — blocks have length L. This notion will 
be essential in the subsequent diversity analysis (Sec. IIVI ). 

In the source phase, relay u,. attempts to decode the received 
symbols from sources belonging to the decoding set S{ur)- 
The users that are successfully decoded at relay Ur are added 
to its retrieval set, denoted by Tl{ur), 72.(u,.) C S{ur), with 
cardinality lu^. Next, in the relay phase, relay Ur transmits a 
relay packet, which is a linear transformation of source 
codeword^ originated by the sources from the transmission set 
T{ur) ~ {ui, . . . , of relay Ur, with T{ur) C TZ{ur)- 

If lu^ < Tiu^, then relay Ur does not transmit anything. In Sec. 
IIVI we show that is an important parameter that strongly 
affects the diversity order 

For example, user 3 attempts to decode the messages from 
users 1, 2 and 5, and succeeds in decoding the messages from 
users 1 and 5 from which a linear transformation is computed. 
Hence, 5(3) = {1, 2, 5}, 7^(3) = T{3) = {1, 5}, /s = ?^3 = 
2. Because the channel between a node and the destination 
remains constant during both source and relay phases, a relay 
has no interest in including its own source message in S{ur). 

Using the transmission set for each relay, the GLNC in Eq. 
(|2]i generalizes to 



(5) 



where the matrices H^^ and iJ^'' are of dimension K x L. 
Hence, each transmitted relay codeword r^^. is a linear trans- 
formation of riu^ source codewords. The superscript Ur in iJ"^ 
indicates that the vector is in general not transformed by 
the same matrix for all relays Ur where Ug £ T{ur). The 

^In practice, increasing the SNR value can be achieved by increasing the 
transmission power of a node, so that both the SNR of the node-to-destination 
channels and channels between non-destination nodes increase. 

"^For conciseness, we do not formulate the equation for channels between 
non-destination nodes. 

'Note that relays u are not allowed to consider relay codewords r^^ for 
inclusion in S{u). As a consequence, the right part of -ffoLNC is diagonal in 
Eq. [7] This restriction was not always applied in the hterature (e.g., (2|), but 
it simphfies the theoretical analysis and code design. 
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overall parity-check matrix H is thus expressed as 



H = 



He 

^^GLNC 



(6) 



where He is block diagonal with Hp on its diagonal, repre- 
senting the channel code, and 

Hi ■■■ H^n. 
Hi ■■■ 



Hglnc 
represents the GLNC. 



Hi 
H2 







Hr, 



(7) 



IV. Diversity analysis of JNCC 

Before passing to the actual diversity analysis, we provide 
the well-known formal definition of the diversity order ||39l . 

Definition 1 The diversity order attained by a code C is 
defined as 

a = — lim — . 

7^oo log 7 

In other words, Pgw oc 7"'', where oc denotes proportional to. 

In the proofs of propositions in this paper, we will often 
use the diversity equivalence between a BF channel and a 
block binary erasure channel (block BEC), which was proved 
in Is), in. A block BEC channel is obtained by restricting the 
fading gains in our model to belong to the set {0, oo}, so that 
a point-to-point channel is either erased or perfect. Denoting 
the erasure probability Pr[Q;„__ = 0] by e, a diversity order d is 
achieved if Pew oc e'' for small e lfT4l . A diversity order of d is 
thus achievable if there exists no combination of d — 1 erased 
point-to-point channels leading to a word error. On the other 
hand, a diversity order of d is not achievable if there exists 
at least one combination of d — 1 erased channels leading to 
word error 

In this section, we present the relation between the diversity 
order d and the parameters — 1, . . . ,mr}, as well 

as between d and the choice of {T{ur),Ur = 1,...,™^}- 
This guides the code design and furthermore, the potential, of 
a linear binary JNCC satisfying some conditions, to achieve 
full diversity, can be verified without performing Monte Carlo 
simulations. 

We first prove that the diversity order is a function of only 
the network coding rate i?„ (Sec. lIV-Al i. We then determine in 
Sec. lIV-Bl the relation between the diversity order d and the set 
{ri„__, = 1, . . . , rrir}, for any linear binary JNCC expressed 
as in Eqs. (|6]l and (|7|i. The set {ri„_ , = 1, . . . , rrir} actually 
determines the maximal spatial diversity that can be achieved 
by cooperation, leading to an upper bound on the diversity 
order In Sec. lIV-Cl we propose a lower bound on the diversity 
order in the case that riu^ = n = 2, which depends on all 
transmission sets {T(ur), u,, = 1, ■ • • , m,.}. In Sec. IIV-DI we 
discuss how the diversity order is affected by interuser failures. 
Finally, in Sec. IIV-EI we briefly comment on the diversity 
order in a layered construction, such as the OSI model. 



A. Diversity as a function of the network coding rate 

We denote the maximum achievable diversity order by dmax- 
We will determine dmax in this section and show that it only 
depends on the network coding rate R„ 



Proposition 1 Under ML decoding, the maximum diversity 
order dmax that can be achieved by any linear JNCC is 



1 + TO,. — nig , if nir > 2ms 



(8) 



Proof: See App. |A] ■ 
Note that the maximal diversity order does not depend on 
L. It can actually be reformulated in the following way: 



p l + (l-H„)(m,.+m,) ^ , if TO, < 2to, 

1 + TOr — (TOs + mr)Rn , if > 2ms 



(9) 



which for to, = to^ = m reduces to the maximum diversity 
order for a standard BF channej^ with m blocks and coding 
rate i?„ QS), Gl, USD. 

Hence, the maximum diversity order does not change when 
the point-to-point channel coding rate Rc.p changes. This 
corresponds with our intuition as the parity bits of the point- 
to-point codes only provide redundancy within one block 
forming a point-to-point codeword, hence these parity bits 
cannot combat erasures which affect the complete point-to- 
point codeword. Another consequence is that the maximal 
diversity order of JNCC cannot be larger than in a layered 
approach, with the same network coding rate. 

In the remainder of the paper, full diversity refers to the 
diversity order being equal to the maximal diversity order, 
d = limax, from (Ull. 

B. Space diversity by cooperation 

We denote the word error rate for each source Us by Pew,us , 
which is the fraction of packets where at least 1 of the K 
information bits from source Us is erroneously decoded at 
the destination. Associated to Pew.u,. we define so that 



P.V 



oc 



for large 7. We have that max„ Pe\ 



X^u ^ew,!is- From Def. [T] it follows that 

d = min du ■ 



(10) 



Denote Us G {!,..., to^}, as the number of times that 
source Us is included in the transmission set of a relay: t^^ = 
£ 7~(ur)), where 11(.) is the indicator function, 
which equals one when its argument is true and zero otherwise. 
Some simple measures can be determined: t,nin = iniiiu, tu^ 
and <av = — ^^^^ — —■ We will show that c?„^ depends on tu^ 
and thus, by Eq. (fTOl i, d depends on tmin- We denote 1 + tmm 
by dii, which we call the space diversity order, as it is the 
minimal number of channels that convey a source message to 
the destination. 



''A standard BF channel is a channel with B blocks of length L, where 
each block is affected by an independent fading gain. The maximal achievable 
diversity order on this channel is given by 1 + [-B(l — Rc)\, where Rc is 
the coding rate (B), {24), (3D- 
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Proposition 2 For any linear JNCC, applied in our system 
model, the diversity order d is upper bounded as 

d < dn = 1 + tmin- 

Proof: We use the diversity equivalence between a BF 
channel and block BEC 15], 161. Assume that the channel 
between source Ug and the destination is erased. Source Us 
is included in at most tu^ transmission sets. Assume that all 
channels between the relays, that include source Ug in 
their transmission set, and the destination are also erased. Then 
the destination does not receive any information on source Ug 
so that it can never retrieve its message. The probability of 
occurrence of this event is e^+'^s, so that Pew,us ^ 6^+*"% 
hence du^ < 1 +tu^- Using Eq. (fTOl i. we obtain Prop. [2] ■ 
Note that the proof of Prop. |2] is based on the assumption that 
relay Ur only considers packets transmitted in the source phase 
for inclusion in S{ur). In the case that relay Ur computes 
its relay packet also based on packets transmitted by other 
relays during the relay phase, the diversity order becomes more 
difficult to analyse. 

In Cor. [T] we propose the conditions on in,in so that the 
space diversity order du is not smaller than the maximum 
achievable diversity order. 

Corollary 1 For any linear JNCC, applied in our system 
model, full diversity can be achieved only if t„m > <?, where 



, if rrir < 2ms 
, if rrir > 2ms 



Proof: The proof follows directly from Props. [T] and |2] 

■ 

Given a GLNC, and thus a choice of T{ur), one can 
verify through Cor. [1] whether full diversity can be achieved. 
However, to get more insight for the code design, we consider 
the simplest case of a network code where the cardinality of 
the transmission set is constant (riu^ = n). 

Corollary 2 For any linear JNCC, applied in our system 
model, with constant ri,„^ — n, full diversity can be achieved 
only if 



n > 
n > 



7n 



n > 



, if mr — ms — m 
if 2ms > "Tir > nT-s 
, if mr > 2ms 



(11) 



Proof: It always holds that tmin < L^avJ and if riu^ ~ n, 
then = From Cor. [T] full diversity can be achieved 



only if 



> q. Because 



> 



we have the 

necessary condition that n > q—. As n is an integer, this 



'lis ' 



> 



Filling in q 



bound can be tightened, yielding n 
from Cor. [1] yields Cor. [2] 

Table J] illustrates Cor |2] showing the set of networks in 
which a certain parameter n is diversity-optimal, which means 
that the choice of n does not prevent the code to achieve 
full diversity. In Sec. |V] we propose a JNCC for n = 2, 
where taking n = 2 is diversity-optimal in all networks 
corresponding to bold elements in Table |I] 
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TABLE I 

Minimal value n for a JNCC with constant = n to maintain 

ITS CAPABILITY TO ACHIEVE FULL DIVERSITY. 



C. A lower bound based on {T{ur)} for n^^ ~ 2 

A certain relay does not help one source only, but a com- 
bination of sources, expressed by the transmission set T{ur) 
for each relay Ur- In this section we provide a lower bound 
on the diversity order, based on the choice of {T{ur),Ur = 
1, . . . , mr}. If this lower bound and the upper bound in the 
previous section are tight, the exact diversity order of JNCCs 
can so be determined, as will be illustrated in Sec. |V] 

Based on T(ur), rng and m,., we construct the {ms + mr) x 
TOs coding matrix AI, where 



M 



1 for Us = I, ■■■ ,ms 

+m,,«, = 1 if Us e 7"(ltr),V 

otherwise 



(12) 



The matrix M expresses the presence of a source-codeword 
in each transmission, i.e., A/i.u^ = 1 if is considered in 
transmission i (i = 1, . . . , and i = ms + 1, • . • , ms + 
mr correspond to the source and relay transmission phases, 
respectively). Therefore, the upper part of M is an identity 
matrix as each source Us transmits its own codeword s^^ in 
the source phase. The matrix M represents what is often called 
the "coding header" or "the global coding coefficients" in the 
network coding literature (see e.g. JT]). 

Consider a block BEC channel where e of the blocks 
have been erased. The indices of the fading gains corre- 
sponding to the erased blocks are collected in the set £ = 
{£i, . . . ,£e},£i S {1, . . . , TTir}). Based on £, we construct 
Ms which corresponds to the subset of transmissions that are 
not erased, i.e., all rows £i (if £i < ms) and ms + £i, for 
i = 1, . . . , e, in M are dropped. We denote the rank of Mg as 
rMe- The set A4{e) collects all possible matrices Mg which 
can be constructed from M if \£\ = e. 

Consider an example for ms = = 3. Assume that 
r(l) = {2,3}, r(2) = {1,3} and T{3) - {1,2}, so that 



M 



1 











1 











1 





1 


1 


1 





1 


1 


1 






(13) 



Next, assume that £ = {1}. Hence, the channel between user 
1 and the destination is erased, so that rows 1 and 4 from M 
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are dropped: 

"0 1 



1 1 

and 3. It can be verified that all matrices G A4{1) 

have rank tmi; = 3. However, there exist matrices Mg G 
M{2) having rank < 3. 

We can now define a metric that depends on {T{ur)}- 



substitution, mentioned in App. |B] corresponds to solving the 
K unknown information bits of s„ via the set of K equations 

H^^^u^ i?:;s„^ +H„^y,„^+„^. (17) 

Us^U 

In Sec. lYl we propose a JNCC where the parity bits of point- 
to-point codes do not have a support in HglnC; so that we 
take ( [TtI i instead of ( fT6b as condition for BP decoding in the 
remainder of the paper 



Definition 2 We define cIm = e* + 1, where e* is the maximal 
cardinality of £ such that = Trig for each g A4{e). 

A simple computer program can compute dM, given 
T{ur),ms and m^- 

Lemma 1 In a JNCC following the form ofEq. (|6| with rus = 
nir and constant ri„^ = n = 2, the metric d^i is at most three. 

Proof: If TTis = rrir and n = 2, then the minimum column 
weight of M is smaller than or equal to three. Erasing the three 
rows where = 1, for a certain Us corresponding to the 

minimum column weight, leads to Mg having at least one zero 
column, and thus r^^ < mg. By Def. |2] dM < 4. ■ 
In the next proposition, we provide a lower bound on the 
diversity order under ML decoding or Belief Propagation (BP) 
decoding 1321 . We denote 



Hp 

Hit'- 



T~Lu,, — 



Hp 



(15) 



which are square matrices of dimension L. 



Proposition 3 Using ML decoding, the diversity order of a 
JNCC following the form of Eq. with constant riu^ = n = 
2, is lower bounded as 

d > dM, 

if the matrices H"^, Us G T{ur), Ur £ {1, . . . , rUs}, have full 
rank. 

Using BP-decoding, the diversity order of a JNCC following 
the form of Eq. (|6| with constant = n = 2, is lower 
bounded as 

d > dM, 

if for each Ur, the set of L equations 



(16) 



Us GT{Ur 



can be solved with BP in the case of only one unknown source- 
codeword vector 

Proof: See App. E ■ 
We can simplify the condition for BP decoding, stated in 
Prop. [3] when we assume that the parity bits of point-to-point 
codes do not have a support in Hglnc- or said differently, 
when the L — K right most columns of the matrices and 
Hy^ are zeroes. In that case, one iteration in the backward 



D. Diversity order with interuser failures 

It is often easier to prove that a particular diversity order is 
achieved assuming perfect interuser channels (see for example 
in Sec.|Vl). Here, we discuss how this diversity order is affected 
by interuser failures. 



Lemma 2 In the case of non-reciprocal interuser channels, 
any JNCC achieves the same diversity order with or without 
interuser channel failures. 

Proof: See Appendix ICl ■ 
In the case of reciprocal interuser channels, the achieved 
diversity order with interuser failures depends on the transmis- 
sion sets {T{ur), Mr = 1, • • ■ 7 rrir}. We propose an algorithm 
to construct {T{ur)} in Sec. W\ and we will then discuss the 
diversity order with reciprocal interuser channels. 



E. Diversity order in a layered construction 

In a layered construction, such as the standard OSI model, 
the destination first attempts to decode the point-to-point trans- 
missions. If it can not successfully retrieve the transmitted 
point-to-point codeword for a particular node-to-destination 
channel, then it declares a block erasure, where a block refers 
to one point-to-point codeword. Denoting this block erasure 
probability by e, we have that e « i 1391 . If for example e 
blocks of length L are erased, then the decoding corresponds 
to solving a set of equations with eL unknowns. 

Standard linear network coding consists of taking linear 
combinations of several source packets. In general, non-binary 
coefficients are used in the linear combinations. Hence, pack- 
ets are treated symbol-wise, which is shown to be capacity 
achieving for the layered construction [|23 ,|. Hence, in Eq. (|5]l, 
the matrices H^^ and H!^'' are replaced by identity matrices, 
which are multiplied with a non-binary coefficient in general. 
A consequence of this symbol-wise treatment is that the 
effective block length of the network code reduces to ras+rrir 
and the set of equations is expressed by the coding matrix 
A/g. At this block length, ML decoding (which is equivalent to 
Gaussian elimination at the network layer) has low complexity. 
Therefore, non-random linear network codes being maximum 
distance separable (MDS) achieve the diversity order dM (Def. 
|2]i. Also note that random linear network codes are MDS codes 
with high probability for a sufficiently large field size l2TI . 
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V. Practical JNCC for riu^ = 2 

In the literature, a detailed diversity analysis is most often 
lacking. Codes were proposed and corresponding numerical 
results suggested that a certain diversity order was achieved on 
a specific network. It is sometimes not clear why this diversity 
order is achieved, and how it would vary if the network or 
some parameters change. In the previous section, we made a 
detailed diversity analysis of a JNCC following the form of Eq. 
(|6]l. However, the utility of for example Prop. [3] is limited to 
JNCCs following the form of Eq. (|6]l with a constant ri^^ = 2, 
which suggests that it is very hard to rigorously prove diversity 
claims in general. However, the modest analysis made in Sec. 
ITV] can be applied in some cases and we will show its utility 
through an example. 

We consider networks with = m.^ = m > 4 and a 
JNCC following the form of Eq. (|6]l with = n = 2 for 
Ur = 1, ■ ■ ■ ,m. We will rigorously prove that a diversity order 
of three is achieved, using the Props, of Sec.|IV] From Table |I] 
it can be seen that this JNCC is diversity-optimal for m = 4 
and m = 5. In Sec. IVIII we provide numerical results for 
m — 5. 

From Table |I] it is clear that restricting n to two is not 
diversity-optimal in larger networks. However, it also has 
some advantages. If n = 2, then every relay just needs to 
decode 2 users, and encoding is restricted to taking a linear 
transformation of only two source packets. Furthermore, taking 
n ~ 2 does not impose infeasible constraints on the number 
of sources in the vicinity of a relay in the case that spatial 
neighbourhoods are taken into account. Next, the theoretical 
analysis is simpler in the case n — 2. Finally, taking n = 2 
allows to reuse strong codes designed for the multiple access 
relay channel, e.g. in ifTOl . ifTTI . 

Besides the diversity order, we indicated in Sec. U that 
scalability is also very important. The JNCC proposed here 
is scalable to any large network without requiring a redesign 
of the code. This means that we provide an on the fly con- 
struction method. The latter is particularly important for self 
regulating networks. As a node adds itself to the network, it 
can seamlessly integrate to the network. Together with the new 
symbols sent by the new node, a new JNCC code is formed 
which still possesses all desirable properties. Finally, note that 
due to the large block length of JNCC, ML decoding is too 
complex and low-complexity techniques, such as BP decoding, 
must be used. 

Hence, two properties are claimed: scalability to large net- 
works and a diversity order of three (which is full diversity in 
some cases) under BP decoding. The JNCC code is presented 
in two steps. First, we present the design of {T(wr)} and 
thus the coding matrix M. In a second step (Eq. (|20]i), we 
specify the matrices H^^ and iJ"'' and we will prove that the 
scalability and the diversity order of three are achieved. 

A. First step: design ofT{ur) 

The transmission sets have a large impact on the 

diversity order For example, in 1(121 . a random construction 
was studied (each relay chooses n = 2 sources at random) 
and it was shown that ]E[t.uJ = 2, but Var[t.uJ = 2 as well. 



so that most probably tniin < 2 and < 3 (Prop. |2]i. So we 
need a more intelligent construction. 

We present an algorithm to determine {T(w,.)}, given 
and rrir, and we subsequently determine the corresponding 
metrics imin and d^. We define the function fm,{x) = 
((.T — 1) mod JTig) + 1 which adapts the modulo operation 
to the range 1 < fmA^) — ^s- 

Algorithm 1 Choose transmission set T{ur) 
for each relay u,. = 1 ^ rrir do 

Set Ui = fniA^r + 1) 
Set U2 = ,f m, [Ur + 2) 

Set S{ur) = {ui, U2} 
if {ui,u2} C TZ{ur) then 

T{Ur) = {ltl,M2} 

else 

T{ur) = {}: relay stays silent 
end if 
end for 

The transmission set T{ur) is expressed via the bottom part 
of M. An example of such a matrix M is given in Eq. (fTSl l 

for ms — nir — 5. 



M = 



(18) 



If a node is added as a source node, it adopts the largest 
source index, TOs + I, and relay-only nodes, with indices larger 
than or equal to + 1, increment their index by one. The 
function fms{x) is updated to the new nis- Note that the 
algorithm corresponds to a deterministic cooperation strategy, 
which avoids extra signalling to the destination regarding the 
code design. 

We first consider the case of perfect interuser channels and 
prove that Alg. [T] yields d = 3 (Cor |3]l. We then consider 
interuser failures and prove that the diversity order is not 
affected (Lemma |3]l. 

Corollary 3 Having perfect links from sources to relays, the 
diversity order of a JNCC, with rris = and with trans- 
mission set constructed via Algorithm [7] achieves a diversity 
order d ~ 'S using BP-decoding, if, for each u^, Eqs. M7^ can 
be solved with BP in the case of only one unknown source- 
codeword vector. 

Proof Because the links between sources and relays are 
perfect, the relays will never stay silent. In the case that m,. = 
rus and = 2, we have that tmin = iav = 2 and so dfl = 3. 

Next, we show that cZm = 3 (and thus, according to Lemma 
[H dM is maximized if rt = 2). Consider \£\ = 2. Without loss 
of generality, consider that 8 = {1,2}. Consider the set of 
equations Af^z = c. Variables 23, ... , z„i^ can be recovered 
via the top m,, — 2 rows of Mg. The two relays ui and U2 
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having source Us in their transmission set (T{ui) and T{u2), 
respectively) are 



of the MARC-JNCC is given by Eq. (A.7) in fTOl 



Ml 



JmA^s - 1),W2 = /m,(Ws - 2). 



Hence, source 1 is included in T{m—1) and T{m), and source 
2 is included in T{m) and Hence, relay transmission 

m — 1 can be used to retrieve source 1 and relay transmission 
m can be used to retrieve source 2, as long as m > 4. Hence, 
A/g has full rank. The generalization to any set £ satisfying 
f I = 2, is straightforward. Therefore, we have that (Im = 3. 

As (in = dM = 3, the proof follows immediately from 
Props. 11] and |3] ■ 

Next, it can be proved that a JNCC applied in our system 
model has a diversity order of three, if it has a diversity order 
of three when all interuser channels are perfect. This is proved 
in general for non-reciprocal interuser channels in Lemma |2] 
and here, we consider reciprocal interuser channels. 



Lemma 3 A JNCC, with transmission set constructed via 
Algorithm\l\ achieves the same diversity order with or without 
interuser channel failures when > 4 or when TOs = m,. = 
TO < 4. 

Proof: See Appendix ID] ■ 
For conciseness, we do not consider the other cases, to^ > 
TO,, < 4. 



B. Second step: JNCC of LDPC-type 

In the first step, we specified {T{ur)} and proved that dn = 
dM = 3 if nir = TOs = TO > 3. According to Cor. [3] a 
diversity order of three is achieved under BP decoding if, for 
each Ur, Eqs. ( fTTI i can be solved with BP in the case of only 
one unknown source-codeword vector. In the second step, we 
specify the sub matrices Hu^, H^^, Wur,Us, to satisfy this 
condition, given that {T{ur)} is constructed according to Alg. 

m 

A simple solution is to replace the K left most columns 
in all K X L sub matrices Hu^, Hu'', Vm^jMs, by identity 
matrices. In this case, the joint network channel coding essen- 
tially reduces to a layered solution: the source-codewords are 
decoded at the relays and simply added according to Eq. (|5]i. 
If the network code is used at the physical layer, it has to deal 
with noise and a more advanced code might be required. 

In the literature, a full-diversity close-to-outage performing 
JNCC for the Multiple Access Relay Channel (MARC) has 
been proposed ifTOl . ifTTI . which is a code in the form of Eq. 
([Til. These codes are such that the set of Eqs. 

Hlsi + H^S2 + Hin = 

can be decoded via BP if only one coding vector Si, S2 or ri 
is erased and the other coding vectors are perfectly known. We 
denote this JNCC by MARC-JNCC. The matrix iJcLNC, marc 



H, 



GLNC, MARC 



lii 2i 
I H 
I 



2i2 
I 

R2 



ri 



(19) 



where Sj = [lij 2ij pj] is the codeword from source j, with 
[lij 2ij] and pj denoting the information bits and the parity 
bits, respectively (j = 1,2); lij and 2ij each contain 
information bits. However, the parity bits pj are not involved 
in i/cLNC, MARC- The matrices Ri, with i = 1,2, 3, are random 
matrices, chosen according to the required degree distributions 
of the LDPC code. To facilitate future notation, we denote 



Hi 



Ho = 



Ri 
I 

I 

-R2 



h: 



Ri 
I 

I 

R2 



si 


S2 


S3 


S4 


S5 


ri 


r2 


ra 


r4 





Hi 


Ho 








H3 

















H'l 


H2 








H3 

















Hi 


H2 








H3 





H2 











H'l 











H3 


Hi 


H^ 
























and Hs ~ R^, so that Hglnc = [^^i H2 H3], where 
Hi = Hi or H'- (it will become clear hereunder which one 
has to be chosen at each relay). In Hi and H2, the first two 
block columns each consist of K/2 columns (corresponding 
to information bits) and the last block column consists of 
L — K columns (corresponding to parity bits from the point-to- 
point codes). The zero block columns indicate that parity bits 
from point-to-point codes have no support in these matrices. 
Now replace all sub matrices Hu^, H,"'' by these matrices, for 
each relay Ur, so that in each block column corresponding to 
information bits, we have a random matrix Ri; this is required 
to conform any preferred degree distribution of the LDPC 
code. For example, iJcLNC can be given by 

ra r4 rs 



ffoLNC =0 0' h'i H2 Q H3 



Hs 

(201 

Each set of rows and each set of columns in H will have 
at least one random matrix, so that any LDPC code degree 
distribution can be conformed. We denote this JNCC by the 
SMARC-JNCC, where S stands for scalable. 

Proposition 4 In a network following the system model pro- 
posed in Sec. \ni\ and using BP, the SMARC-JNCC achieves a 
diversity order d = 3. 

Proof Consider the set of K equations 

HsVur = HiSui +i?2S„2, {ui,U2} e T{Ur). (21) 

In ifTOl . it is proved that this set of K equations can be solved 
using the matrices proposed above. We provide another more 
simple proof here. Consider a block BEC. Because Hi and H2 
are upper- or lower triangular, with ones on the diagonal, the 
unknown K information bits can be retrieved using backward 
substition, hence it can be retrieved with BP as well. 

By Cor [5] and Lemma [3] the SMARC-JNCC achieves a 
diversity order d ~ 3. ■ 

'The attentive reader will notice that the first two block rows in Eq. (A.7) 
in 1101 are not used here. These block rows are only necessary if a source is 
helped by one relay only and no point-to-point codes are available, which is 
not the case here. 
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Note that the information bits of a source need to be split 
in two parts: bits of the type li and 2i. This allows the 
introduction of the matrices i?i and R2 in Eq. [19] so that all 
information bits have a random matrix in their corresponding 
block column in the parity-check matrix. Now, the LDPC code 
can conform any degree distribution. 

VI. Lower bound for the WER 

To assess the performance of the SMARC-JNCC we need 
to compare it with the outage probability limit (Sec. IVI-Al i. 
We show that the outage probability limit is not always tight 
and we propose a tighter lower bound, which is presented in 
Sec lyTBl 

A. Calculation of the outage probability 

The outage probability limit is the probability that the 
instantaneous mutual information between the sources and 
sinks of the network is less than the transmitted rate. The 
outage probability is an achievable (using a random codebook) 
lower bound of the average WER of coded systems in the limit 
of large block length Bl. |[T3l. l33l. 

For a multi-user environment, two types of mutual infor- 
mation are considered. First, it is verified whether the sum- 
rate, Rc in this case, is smaller than the instantaneous mutual 
information between all the sources and the sink. Then, it 
is verified whether each individual source rate, — in this 
case, is smaller than the instantaneous mutual information 
between the nodes, transmitting information for this source, 
and the destination. The outage probability for the MARC was 
determined in ifTOl . ll20l using the method described above. 

The outage probability is 

-Pout = P(^^out)j 

where fout is denoted as an outage event. Similarly as in ifTOI . 
an outage event is given by 

E™ =1 1{Su. ; D) + B^J{R^^ ; D 



Rc > 



rUs + rrir 



rUs ~ nis + rrir 



where 



The terms I {Si, D), I{Ri; D) and I [Si] Rj) are the instanta- 
neous mutual informations of the corresponding point-to-point 
channels with input x G { — 1, 1}, received signal y = aiX + w 
with w ^ CAf{0, ^), conditioned on the channel realization 
ai, which are determined by applying the formula for mutual 
information ISll. BOl: 

I{X; Y\a^) = 1 - EY\{x=i,a^} {logs (1 + cxp [-iya.j])} , 

where ^Y\{x=i,ai} is the mathematical expectation over Y 
given X ^ 1 and a,;. 

We now consider the outage probability of a layered con- 
struction, such as the standard OSI model, where the destina- 
tion first decodes the point-to-point transmissions, declaring a 
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13 
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Fig. 1. The outage probabilities of JNCC and a layered construction are 
compared. The spectral efficiency is /? = 3/7 bpcu. 



block erasure if decoding is not successful. For the network 
code, we assume a maximum distance separable (MDS) code, 
which is outage-achieving over the (noiseless) block-erasure 
channel 1141 . That is, any riig correctly received packets suffice 
for decoding. Accordingly, an outage event for the layered 
construction, denoted as Eout.i is given by 



Ur. = l 



where 



and 



1 - Es,^^ + I] (1 - Er,j) = 



Er,u^ = 1 - BuMHRj-.D) > Rc.p] 



The outage probability for JNCC and a layered construction 
are compared in Fig. [T] for nis = m,. = 5, coding matri5i0 
M given in Eq. ( fTSl l and Rc,p = 6/7. The overall spectral 
efficiency is R = 3/7 bpcu, so that Eb/No ~ -y. The 
main conclusion is that the difference between both outage 
probabilities is only IdB. Hence, on a fundamental level, the 
achievable coding gain by JNCC with respect to a standard 
layered construction is small for the adopted system model. 

B. Calculation of a tighter lower bound on WER 

According to information theory, the outage probability is 
achievable, where the proof relies on using random codebooks. 
However, the nature of the JNCC protocol largely deviates 
from a random code. For example, the parity bits correspond- 
ing to the point-to-point codes are forced in a block diagonal 
structure in He (see Eq. |6]l, which is not taken into account in 
the outage probability limit. In fact, in Prop. [1] it was proved 
that the maximal diversity order does not depend on Rc but on 

*The coding matrix expresses the transmission sets for each relay, which 
is required to determine the outage probability. 
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CND 



H, 



CND 




CND 



AGLNc(a;) 



Fig. 2. The depicted part of the factor graph (using a Tanner notation) 
illustrates that a bit node (bit i on the figure) is essentially connected to two 
sets of check nodes, corresponding with He and -ffoLNC. respectively. A set 
of check nodes is denoted as CND for check node decoder. The LLR-value 
coming from the CND corresponding with He is denoted as Lc- The LLR- 
value corresponding with the channel observation is denoted as Lots- 



Rn, which is not taken into account in the outage probability 
limit. Therefore, we argue that the outage probability limit is 
in general not achievable by a JNCC, which is illustrated by 
means of an example. 

Consider a network with = = 3. The adopted 
point-to-point codes have coding rate i?c,p = 0.5, so that 
Rc = 0.25. We take riu = 2 and adopt the coding matrix M, 
given in Eq. (fTlT l. Because of the small coding rate Rc, the 
outage probability achieves a diversity order of three (Fig. |4|i. 
However, it follows from Prop. [T] that dniax = 2. We therefore 
propose a new lower bound, which takes into account the 
point-to-point codes. 

A bit node is essentially protected by two codes: a point- 
to-point code (He) and a network code (i/cLNc), which is 
illustrated on the factor graph |26!| representation (a Tanner 
notation ll38l is adoptedjj of the decoder (Fig. |2]l. Usually, 
both codes are characterized by separate degree distributions, 
denoted as {Xc{x) , pc{x)) and {Xglnc{x), Pglnc{x)) for 
and iJcLNC, respectively. 

The new lower bound assumes a concatenated decoding 
scheme. At the destination, first the point-to-point codes are 
decoded and then soft information is passed to the network 
decoder This is illustrated in Fig. [3] where the soft information 
is denoted by the log-likelihood ratio (LLR) £obs'- Note that 
the bit node of bit i is duplicated to be able to clearly indicate 
^obs' ■ Applying the sum-product algorithm (SPA) on this factor 
graph or the original factor graph (without node duplication) 
is equivalent (see BTl for a background on factor graphs and 
the SPA). The LLR L^bs' can be viewed as a new channel 
observation as it remains fixed during the iterative decoding 
of the network code (-ffcLNc)- The maximum rate that can be 
achieved by the network code is given by 



1 



; Lobs') + ^ Bu-^I[Ru/tL. 



The terms I{Si]Lo\j^') and I{Ri; Law) are the mutual in- 
formations between the channel input x E { — 1,1} and the 

'For a specific instance, the parity-check matiix can be graphically 
represented by a bipartite graph, denoted as a Tanner graph. The graphical 
Tanner graph representation is equivalent to the factor graph, which can be 
used for decoding. 




CND 




Fig. 3. The bit node in Fig. |2] can be duplicated with a single edge between 
both nodes as shown in this figure. The LLR Lobs' is sum of all incoming 
LLR-values from the left, and contains the soft information which is passed 
to the network code decoder in a concatenated coding scheme. 



random variable Lobs' . conditioned on the channel realization 
ai, determined by applying the formula for mutual information 
0, iOl, i.e., /(X; Lobs' is 

The density of the random variable Lobs' can be obtained by 
means of density evolution |[35l , given the degree distributions 
of the point-to-point code, or by means of Monte Carlo 
simulations, given the actual factor graph of the point-to- 
point code. Both approaches yield to the same results in our 
simulations. 

Similarly to the conventional case, an outage event, denoted 
as fout,2 is given by 



^-'^L„,,-\{x = l,a,} \ log2 1 



i'out,2 — 



Rn > 
I{S 



IZZU /(s^.,;i„b.')+E^;=l s„,/(fl„,,;L„b/) 



> 



rrts -{-rrij. 



Note that the network coding rate is used instead of the overall 
rate Rc, which corresponds to Prop. [U 

The tight lower bound presented here is a valid lower bound 
if the point-to-point codes are first decoded, followed by the 
network code, without iterating back to the point-to-point 
codes. 

Let us now go back to the small network example with 
nis = nir = 3, considered in the beginning of this section. Fig. 
|4] compares the conventional outage probability (Sec. IVI-AI ) 
with the tighter lower bound proposed here. As mentioned 
before, the conventional outage probability has a larger di- 
versity order than what is achievable, while the tighter lower 
bound only achieves a diversity order of two. We are seeing 
a 3dB difference at an outage probability of 10^''. To assess 
the performance of the network code only, given a certain 
point-to-point code, the WER of the SMARC-JNCC should 
be compared with the tight lower bound presented here. In the 
subsequent sections, we always include both lower bounds. 

VII. Numerical results 

In this section, we provide numerical results for the 
SMARC-JNCC. We will clarify the proposed techniques 
on an illustrating network example, where m.g = rrij. = 5 
(Fig. |5]). We use the same network example as in ||2l, lfT2l 
so that a comparison is possible. For simplicity, we assume 
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Fig. 5. The network example that will be used in this document is illustrated. 
The solid lines represent interuser channels, the dashed line is the channel to 
the destination. Only the channels from the perspective of user 1 are shown 
for clarity, but all other users see equivalent channels. 



non-reciprocal interuser channel in the simulation results. 
Note that in the case that nis > 4 and Alg. [T] is used to 
construct {T{ur),Ur = l,...,m,.}, reciprocity is irrelevant 
for our proposed code, as it applies that i ^ T{j) if j G T(«). 

We compare the error rate performance of the SMARC- 
JNCC with the outage probability limit and the tighter lower 
bound, which are presented in Sec. I VII and with standard net- 
work coding techniques (using identity matrices in i^cLNc) and 
a layered network construction (also using identity matrices in 
Hglnc, and where, at the destination, the network code is only 
decoded after decoding all point-to-point codewords separately 
and taking a hard decision). 

The point-to-point code used in the simulations is an irregu- 
lar LDPC code ||35]| characterized by the standard polynomials 
X{x) and p{x) 



2 



dc 

E 

i=2 



PiX 



where \{x) and p{x) are the left and right degree distributions 
from an edge perspective. The coefficients Ai and pi are the 
fraction of edges connected to a bit node and check node, 
respectively, of degree i. The adopted point-to-point code is 
fetched from lfT6l . has coding rate i?c,p = 6/7 and conforms 
the following degree distributions: 



Aa = 0.173, A3 = 0.223, A4 

P24 = 0.96, P25 = 0.04. 



0.095, As = 0.51 



A. Perfect source-relay links 

We start by assessing the performance of i/cLNC, the bottom 
part of Eq. ( |20l i, which determines the diversity order. There- 
fore, we assume perfect links between sources and relays. 
Hence, the channel model is the same as described in Sec. 
Uni with the exception of the interuser channels, which are 
assumed to be perfect (no fading and no noise). The parameters 
used for the simulation are K ~ L ^ 900, nig = irir = 5 (so 
that N = \QK = 9000), where N is the block length of the 
overall codeword. The overall spectral efficiency is R = 0.5 
bpcu, so that E^/Nq = 27. 

Fig. |6] shows that a diversity order of 3 is achieved for 
SMARC-JNCC, which corroborates Cor. [3] It performs at 
2.5dB from the outage probability (because no point-to-point 
codes are considered, only the conventional outage probability 
can be calculated), which may be improved by optimizing the 
degree distributions. We also show a JNCC, where all subma- 
trices Hu^, H^J , yur,Us are replaced by identity matrices, 
denoted as the I- JNCC. Finally, we show an I- JNCC with 
irregular with coding matrix AI, given by 



M = 



(22) 



It is clear that, even without optimizing the SMARC-JNCC, 
there is a benefit in terms of coding gain compared to the 
I-JNCC. 

B. Rayleigh faded source-relay links 

Now, we assess the performance of the complete parity- 
check matrix H of the SMARC-JNCC. We use the channel 
model as described in Sec. |III] Hence, all links have the 
same statistical model and the average SNR is the same as 
for all channels. The parameters used for the simulation are 
K = 606, Rc,p = 6/7, L = 707, m, ^ rur ^ 5 (so 
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Fig. 7. The word error rate of the SMARC-JNCC is compared to that 
of the I- JNCC and a layered construction, assuming Rayleigh faded source- 
relay channels. The reference curve is the performance of the SMARC-JNCC 
assuming perfect source-relay channels (Sec. IVII-A) . 
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Fig. 8. The word error rate of the SMARC-JNCC is compared to that of the 
1-JNCC. assuming Gaussian source-relay channels. The reference curve is the 
performance of the SMARC-JNCC assuming perfect source-relay channels 
(Sec. IVll-At . 



that N = lOL = 7070). The overall spectral efficiency 
is i? = 3/7 bpcu, so that Ei,/Nq ~ Because the 

simulation time would be very large if every point-to-point 
source-relay link has to be decoded separately, we made 
an approximation. The word error rate of the point-to-point 
code when transmitted on a channel with fading gain a is 
smaller than 10""' when a^7 ~ 5.5dB. Therefore, we assumed 
that a relay had correctly decoded the source-codeword if 
a^j > 5.5dB and not otherwise. We also add the performance 
of the SMARC-JNCC from Sec. IVII-AI coiTesponding to 
perfect source-relay links and R = 0.5 bpcu, as a reference 
curve (note that the reference curve corresponds to a larger 
spectral efficiency - the coding rate Rc is larger - than for the 
other curves, which slightly disadvantages the reference curve 
in terms of error performance). 

Fig. |7] shows that a diversity order of 3 is still achieved, 
which corroborates Prop. H] In addition, two main conclusions 
can be made. First of all, the loss due to interuser failures is 
6.5dB, which is very large. Secondly, the benefit in terms of 
coding gain of the SMARC-JNCC compared to the I-JNCC 
is considerably decreased, compared to Sec. IVII-AI which 
corresponds to the small horizontal SNR-gap between the 
outage probabilities of a layered and joint construction. Also 
note that the tighter lower bound using density evolution, is 
close to the conventional lower bound in this case. Finally, the 
WER performance of a layered construction is shown, which 
coincides with that of the I-JNCC. 

C. Gaussian source-relay links 

We test again the complete parity-check matrix H of the 
SMARC-JNCC, now assuming that the source-relay links are 
Gaussian, having additive white Gaussian noise only, without 
fading; fading occurs on the source-destination and relay- 
destination links only. We assume that the average SNR is the 
same for all channels. The parameters used for the simulation 
are the same as in Sec. IVII-BI 

Fig.[8]shows that in the case of Gaussian interuser channels, 
the loss compared to perfect interuser channels is very small. 



Furthermore, the performance of the I-JNCC has improved a 
lot in comparison with Sec. IVII-AI where i?GLNC only was 
used. The degree distributions causing the poor coding gain 
of the I-JNCC in Sec. IVII-AI are considerably improved by 
the point-to-point codes. 

VIII. Conclusion 

We put forward a general form of joint network-channel 
codes (JNCCs) for a wireless communication network where 
sources also act as relay. The influence of important parameters 
of the JNCC on the diversity order is studied and an upper and 
lower bound on the diversity order are proposed. The lower 
bound is only valid for the case where the number of sources 
is equal to the number of relays, and where each relay only 
helps two sources. 

We then proposed a practical JNCC that is scalable to 
large networks. Using the diversity analysis, we managed 
to rigorously prove its achieved diversity order, which is 
optimal in a well identified set of wireless networks. We 
verified the performance of a regular LDPC code via numer- 
ical simulations, which suggest that as networks grow, it is 
difficult to perform significantly better than a standard layered 
construction. 
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Appendix 

A. Proof of Prop. [7] 

The maximal diversity order can be derived using the di- 
versity equivalence between a block BEC and a BE channel 
ifSl , 161 . Assume a block BEC, so that a block or is 
completely erased or perfectly known. Consider the case that 
ei blocks of length 2L and 62 blocks of length L have been 
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erased, where e = ei +62 is the total number of erasures, ei < 
nis and 62 < nir — rag. Hence, the number of unknown bits 
is equal to ei2i + 62^. Considering the structure of H from 
dill containing the block-diagonal matrix He, it follows that the 
ei2i+e2L erased bits appear in only (2ei+e2)(L— A')+mr/\ 
of the available (m^ + mr)L — nisK parity equations, i.e., 
(2ei + 62) (i — K) equations involving He and all rrirK 
equations involving iJcLNC- Hence, the unknown bits can be 
retrieved only if there are sufficient linearly independent useful 
equations. This yields the necessary condition: 



nir > 2ei + 62- 



(23) 



Denoting by e = ei + 62 the total number of erased blocks, 
the largest value Cmax of e for which ei and 62 satisfy (|23] i for 
all ei < TTis and 62 < rrir — is given by 

[^J rrir < 2ms 

mr — TO, mr > 2to, 



(24) 



Hence, dma, 



1, yielding Prop. [T] 



B. Proof of Prop. \3\ 

Before we present the actual proof, we first propose two 
lemmas. 

Lemma 4 Any binary axb matrix S, a > b, where all rows 
have weight 2 cannot have full rank b. 

Proof: If a matrix has full rank, there is no vector z 7^ 
such that Sz = 0. However, if S has row weight 2, then 
51 = 0, where 1 corresponds to a column vector with each 
entry equal to 1. ■ 
Consider now a column vector of b unknown variables 
z and a set of constraints on these variables, which are 
stacked in S so that Sz = c, where c is a column vector of 
known constants. In general, solving S for z corresponds to 
performing Gaussian elimination of S. However, under some 
conditions, this simplifies to backward substitution. 

Lemma 5 If a binary axb matrix S, a > b, has full rank b 
and maximal row weight of 2, Gaussian elimination simplifies 
to backward substitution. 

Proof: Without loss of generality, we eliminate all redun- 
dant (linearly dependent) rows in S to obtain a square matrix 
of size b. By Lemma 21 there must be at least one row in S 
with unit weight to have full rank. Starting from this known 
variable, we can solve for a further variable in z at each step 
as the row weight is smaller than or equal to 2. 

Assume that this backward substitution procedure cannot 
be continued until all variables are known. That is, after 
successive decoding, there are k rows consisting of a com- 
bination of Zij, + Zjj, where neither Zi^ nor Zj^ are known. 
We split the matrix S into two parts; S'unknown G {0,1}*^^™= 
and S'known G {0, The former comprises the rows 

involving only unknown variables (note that the weight of each 
row of 5unknown IS 2). The latter consists of the rows involving 
only known variables. If the number of unknown variables is 
equal to k, then the rank of ^unknown must be equal to k which 



is impossible by Lemma |4l So, the matrix S was not full rank 
which contradicts our assumption. If the number of unknown 
variables is smaller than k, then there were redundant (linearly 
dependent) rows in S'known which contradicts the assumptions 
again. We conclude that the procedure only fails if 5* does not 
have full rank. ■ 
To prove Prop. |3l we use the diversity equivalence between 
a block BEC and the BF channel. In a block BEC, the channel 
equation (jUi simplifies to 



Us 
Ur 



1, . . . ,TOs 

l,...,mr 



(25) 



where = when the channel is erased and = 1 otherwise. 
Hence, ei = if i G £ and ei = 1 if i E £, where £ is the 
complement of £. 

Source-codewords s; can be retrieved from the transmis- 
sions in the source phase if ei = 0. Decoding the other source- 
codewords at the destination is performed through the parity- 
check matrix H (Eq. We split H in two parts: 



left 



H 



right 



(26) 



have msL and 77i,.L columns. 



where TJieft and iJright 
respectively. We also define s = [s[ 
r = [r^ . . . Y^^y . As i7x = 0, we have that 

i/leftS = -ffrightl". 



-T iT 



and 



(27) 



As we consider a block BEC, some transmissions are perfect. 
As in App. |A1 consider the case that ei blocks of length 
2L and 62 blocks of length L have been erased, where 
e = ei + 62 = \£\ is the total number of erasures, ei < to^ 
and 62 < mr — to^. Considering the structure of H from ^ 
containing the block-diagonal matrix He, it follows that the 
ei2i+62i erased bits appear in only {2ei+e2){L— K)+mrK 
of the available (m^ + mr)L — msK parity equations, i.e., 
(2ei + 62) (i — A') equations involving H^. and all TO,.i^ 
equations involving Hglnc- Next, (ei + e2)K from the TO,.i^ 
equations involving i?GLNC cannot be used to solve erased bits 
in s as these equations always have at least two unknowns. The 
overall set of equations to decode s thus becomes 



or, using the notation from ( fTsl l. 



\f Us E £ 

y use£ (28) 

y Ur G £, 



y Us € £ 



(29) 



where y[ = (BPSK modulation). We can stack the 

coefficients of all elements in s in a matrix Hg. For example, 

ifms=mr^3,£ = {!}, T{2) = {1, 3} and r(3) = {1, 2}, 
then 



Hs 



si 


S2 


S3 





I 











I 


Hi 





Hi 


Hi 


"2 






(30) 



It is now easy to see that Mg, as defined in Sec. 
HV^ is closely related to Hs'. [Afflij = 1 if 
[Hs]{i-i)L+i...iL,u~i)L+i...jL 7^ and [M£]ij = if 
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[Hs]{i~l)L+l...iL.{j-l)L + l...jL = 0- 

If |f I < cJa/ — 1, then Ms has full rank, according to Def. |2] 
As established in Lemma |5] the set of equations represented 
by M£ can be solved using backward substitution. This means 
that at each iteration, there is an equation with only one 
unknown. Consider a particular iteration and denote the index 
of the unknown by u. In Hg, this corresponds to an equation 
with an unknown source-codeword vector Su of the type 







(31) 



1 Hu-s^ 

L Us 

or of the type s„ y^. 

Under ML decoding, we obtain what was claimed if the 
matrices H"'', Ug G T{ur), Ur G {!,..., rrir} have full rank. 
Under BP decoding, we obtain what was claimed if, for each 
Ur, the set of L equations ( [3T] i can be solved with BP in the 
case of only one unknown source-codeword vector s^. 

C. Proof of Lemma \2\ 

A relay may not succeed in successfully decoding the mes- 
sage from a source, denoted as a failure. There are ~m in- 
teruser channels, which all have a probability of failure. Hence, 



m — m 



there exist X^Tio 



different cases, where each 



m — m 
i 

case corresponds to a combination of failures and successes. 
We denote the case where all interuser channels are successful 
as case 1. 

Using Bayes' law, the error rate can be split: 



Pew = P(case i)P(ew|case i). 



(32) 



Defining the diversity order corresponding to each case as 

dc,^ = - lim^_^oo '°sP{c^^<^^)P^{^^\c^^^ '0 ^ it follows that the 
overall diversity order is d = mini dci- 

The probability of / failures on independent interuser chan- 



nels is proportional to 1391 . so that for this case i, 

loffPfcase i) loffPfewlcase i) ,,,, 
4,^ = - lim -^-^ - hm -^-^ — ^ '- (33) 

7-foo log 7 7->oo log 7 



/ — lim 

7— >CJO 



P(ew|case i) 
log 7 



(34) 



The diversity order in the case of perfect interuser channels 
(/ = 0) is del. That is, the error-correcting code can bear 
dc^i — l erasures on node-destination links. Hence, dc,i > (ic,i 
only if P(ew|case i) < df^-f , or, all information can still 
be retrieved at the destination, given that / interuser channels 
and dc,i — / — 1 node-destination channels are erased. Let us 
check whether this is true for all /. 

A relay stays silent if it cannot decode all source codewords 
corresponding to its transmission set. If there are / interuser 
failures, there are at most / relays which stay silent in the 
relay phase. This corresponds to at most / additional node- 
destination erasures adding to the assumed d^.i ^ / ^ 1 
already erased node-destination channels, yielding a total of 
at most dc_i — 1 erased node-destination channels, which can 
be supported by the code, by the definition of dc.i- 



D. Proof of Lemma \3\ 

In the case that nis > 4 and Alg. [T] is used to construct 
{T{ur),Ur ~ 1, . . . , m,.}, reciprocity is irrelevant for our 
proposed code, as it applies that i ^ T{j) if j G 7~(i). Hence, 
if nis > 4, the proof given in App. |C]is always valid. 

Now consider the case that dc^i — 2, which corresponds 
to 77is = rrir ~ m < A (see Prop. [T]). In the case of 
/ = 1 interuser channel, c?c,i is always larger than one, 
because P(ew|case i) < ^ as at least one channel, the source- 
destination channel, needs to fail to loose the corresponding 
information bits. 

Finally, consider the case that nis = "mr = m = A and 
thus del ~ 3. Hence, in the case of no interuser failures, the 
code can support two node-destination failures, corresponding 
to four erased transmissions from two nodes, in the source 
phase and in the relay phase. Reciprocity is relevant as i G 
r(j) if j G Tii) for is (1,3) and (2,4). Because 

P(ew|case i) < we only have to consider the case that 
/ = 1, denoted as case i in general. Hence, in the case 
that the interuser channel between sources one and three or 
two and four have been erased, relays one and three or two 
and four, respectively, stay silent. Note that the transmission 
sets from the remaining active relays are disjoint when Alg. 
[U is used, and because n = 2, they support all sources 
Ms = 1, . . . , 4. If one node-destination channel is consequently 
erased, which corresponds to at most two transmissions, the 
destination has to recover the information bits from the erased 
source-codeword. Because relay Ur cannot have Ur in their 
own transmission set T{ur), the erased relay codeword does 
not contain any information on the erased source-codeword, 
which implies that the information is in the remaining relay 
codeword. Hence, we have that P(ew|case i) < or by ( |34] |. 
dc,i > 3. In other words, interuser failures do not decrease the 
diversity order. 
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